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Abstract 

It has long been appreciated that transport properties can control reaction kinetics. This effect 
can be characterized by the time it takes a diffusing molecule to reach a target - the first-passage 
time (FPT). Although essential to quantify the kinetics of reactions on all time scales, determining 
the FPT distribution was deemed so far intractable. Here, we calculate analytically this FPT 
distribution and show that transport processes as various as regular diffusion, anomalous diffusion, 
diffusion in disordered media and in fractals fall into the same universality classes. Beyond this 
theoretical aspect, this result changes the views on standard reaction kinetics. More precisely, we 
argue that geometry can become a key parameter so far ignored in this context, and introduce the 
concept of " geometry-controlled kinetics" . These findings could help understand the crucial role of 
spatial organization of genes in transcription kinetics, and more generally the impact of geometry 
on diffusion-limited reactions. 



PACS numbers: 



It is known on general grounds that reaction kinetics can be influenced by the transport 
properties of the reactants [TH3]. Actually, the transport step, before reactants meet and 
eventually react, even becomes limiting in the case of confined systems such as cells or cell 
subdomains where a small number of reactants are involved [IH9]. In such systems, the 
first step in estimating the kinetics of reactions consists in evaluating the properties of first- 
encounter between reactants. Quantitatively, this amounts to calculating the distribution of 
the time it takes a diffusing molecule to reach a target site - the first-passage time (FPT) 
distribution. While this quantity is well-known in quasi-one dimensional or unconfined 
geometries (see [10J for a review), determining the FPT distribution appears to be intractable 
in the realistic situation where the diffusing molecule is confined within a finite domain [TT] . 
A first estimate of the effect of geometrical parameters of confinement on this search time 
is given by the mean of the FPT. This has recently been calculated, and a linear scaling 
with the volume has been demonstrated [T2TH5] . However, as soon as several time scales are 
involved, the kinetics can not be determined by the mean of the FPT only, and the entire 
distribution is needed [To^ [T7|. 

Gene transcription provides an extremely important example - which we shall repeatedly 
invoke in what follows - of reactions involving a small copy number of (or even single) reac- 
tants confined within a small (micron-sized) domain, and whose kinetics must be precisely 
regulated to fulfill vital cell functions. Interest in the question of how geometrical parameters 
impact on the kinetics of such transcriptional reactions and could act as regulatory factors 
has recently increased, mainly because of new experimental tools that enable the observation 
of the real-time production of proteins at the single molecule level [I] . These technics, which 
give access to the spatial organization of the genetic material, have revealed strong correla- 
tions between the spatial locations of successively activated genes, both for prokaryotes [18] 
and eukaryotes [19J. Indeed, it has been found that successively activated genes are often 
colocalized, that is located in the very same nuclear regions. These observations raise the 
question of the importance of geometrical parameters in transcription kinetics, which has 
remained so far widely unanswered. 

In the broader context of chemical reactions in confinement, we are interested here in 
the following questions : (i) How does the FPT distribution depend on the volume of the 
confining domain? (ii) How does it depend on the initial position of the diffusing molecule, 
and (iii) is this geometric dependence an important factor, which could potentially control 
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the kinetics ? 

Note that the influence of confinement and crowding effects on biochemical reactions has 
already been studied (see {201 l2"Tj for reviews) on the basis of a thermodynamical treatment 
of reaction kinetics. While this approach is well suited to the case of a large number of 
reactants, it does not provide the dependence of the kinetics on the geometrical parameters 
mentioned above (such as the initial position of the reactants), which involve the individual 
nature of the reactants and their dynamical properties. 

In this work, we calculate analytically the distribution of the FPT at a target T for 
a diffusing particle released at a starting point S (see figure [T]) and quantitatively answer 
questions (i)-(iii) above. We highlight universal laws of the FPT distribution as function of 
the volume TV of the confining domain and of the distance ST = r and show that two regimes 
emerge. More precisely, we find that the key criterion is the compact vs non compact nature 
of the diffusion process, to be defined mathematically below. In the non compact case, which 
physically corresponds to a diffusing molecule which "sparsely" explores its environment and 
leaves unvisited regions (typically a molecule diffusing in a dilute solution), we show that 
the kinetics is widely independent of the starting point. On the opposite compact case of 
a diffusing molecule which "densely" explores its environment (for example a molecule in a 
very crowded medium), the position of the starting point strongly influences the search time 
of the target, which leads us to introduce the concept of "geometry-controlled kinetics". In 
the context of gene transcription, this result implies that the kinetics of activation of a gene 
T by a transcription factor (TF) can be orders of magnitudes faster if the TF is released 
from a site which is colocalized with (i.e. in the vicinity of) T, as compared to the case 
where the TF is released from a remote site. 

Results 

We consider a Markovian random walker of position r(t), whose dynamics is characterized 
by the dimension of the walk d w , defined by the scaling of the mean squared displacement 
(r 2 (i)) oc t 2 l dw . The walker is confined in a domain of N with reflecting walls. Additionally, 
we assume that the medium is of fractal type, so that its characteristic linear size R scales 
as R oc N 1 / d f , where df is the fractal dimension [22]. We are interested in the distribution 
P(Tts) of the time it takes a walker starting from the site S to reach for the first time the 
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target site T located at a distance r from S. 

We start from the backward equation satisfied by the probability P(T TS = t) in discrete 
time t PU 

P(T TS = f) = X) v>ssP{T TS =t-l), (1) 

j 

obtained by partitioning over the first step of the walk, where it?y stands for the transition 
probability from site j to site i. It is shown in Supplementary Information (SI) that this 
equation, after Laplace transform, leads to the following hierarchy satisfied by the moments 
(T Tj ) of the FPT : 

< T ts> = (^)(-l) fc+ M(#Tr-#TS^ (2) 

Here 

oo 

% = EW) - ( 3 ) 

n=0 

and Wyj(n) denotes the propagator of the walk, ie the probability to be at site j at step n 
starting from site i, while Wj tat is the probability to be at site j in the stationary state. 

At this stage, the hierarchy of equations ([2j remains formal, as it involves the unknown 
functions Hji, and does not allow an explicit determination of the FPT distribution. How- 
ever, this difficulty can be circumvented by taking the large volume limit and by considering 
the rescaled time 9 = T T5 /(T) T , where (T) T = 52 3 {T T s)Wg* stands for the mean FPT 
to the target site r^, averaged over the initial position. Note that here we implicitly assume 
that (T) T (as well as its disorder average in the case of disordered systems to be discussed 
below) is finite. Actually, a detailed analysis of Eq.Q allows us to show in SI that the 
distribution of the rescaled variable 9 takes in the large volume limit the following general 
form 

Gts(0) = (1 - n TS )5(9) + U TS ^(9) (4) 

where the Dirac 5 function corresponds to trajectories hitting the target without reaching 
the boundary within a time of order r dw much smaller than (T) T . The geometrical factor 
1 — Uts can be interpreted as the weight of these trajectories. Similarly, the contribution 
n.Tsi>{0) accounts for trajectories reaching the boundary before the target. Note that the 
dependence on the starting point lies entirely in the geometrical factor Uts, whereas the 
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time dependence is contained in the scaling function ifj. The geometrical factor Uts and the 
rescaled variable 9 are explicitly determined in SI and their scaling with the volume N and 
the distance r are obtained under the standard scale-invariance assumption of the unconfmed 
propagator W°°(n) oc n - d f/ d ^f (| r . — rj\/n l ^ dw ) [22]. Actually, the FPT distribution given 
by Eq. Q falls into a few universality classes, defined according to a purely geometrical 
criterion as detailed below. 

In the case of non compact exploration defined here by d w < df [23], where the mean 
number of distinct sites visited by the walker in absence of confinement grows linearly with 
the number of steps, we find : 



where a is a lattice dependent constant of order 1. Note that the linear dependence on N 
of the scaling variable (T) T is the same as found in [T5] for the mean FPT (Tts), and in 
particular does not depend of the dimensions df and d w . This result includes as a special 
case regular diffusion in 3D given in [24J for which d w = 2 and df = 3. Strikingly, the 
exponential form of ip holds for any dimensions such that d w < df. Note that while the FPT 
distribution is a mere exponential of weight one in the limit of r larger than the step size, it 
departs significantly from this distribution if the starting point is close to the target. 

In the opposite case of compact exploration d w > df [23], where the mean number of 
distinct sites visited by the walker in absence of confinement grows slower than linearly with 
the number of steps, further hypothesis on the unconfmed propagator are needed to estimate 
the relative importance of the terms involved in Eq. ([2]). Making use of the O'Shaughnessy- 
Procaccia operator [25] to evaluate the large volume behaviour of Hij in ([2]), we find that 
the FPT distribution obeys the generic form of Eq. Q with 



where 9 > 0, v = df/d w and «o < «i < ... stand for the zeros of the Bessel function J_ v . 
Strikingly, the scaling with N of the scaling variable 9 is no longer given by the mean FPT 
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(Tts), which clearly indicates that several time scales are involved in the problem. The 
interplay between these time scales leads to a non trivial family of universal non exponential 
scaling functions, parametrized by d w and df. The geometrical factor strongly depends on 
the source-target distance r, and in particular is very small for r <C R, as opposed to the 
non compact situation. We add that the marginal case d w = df, which is compact according 
to the definition given in [23], corresponds to an exponential scaling function ip as given by 
Eq. (jHJ), with logarithmic corrections in the scalings of (T) T and Uts with r and N (see SI). 

Eqs([5])-(|6| fully define universality classes of FPT distributions in confinement. Addi- 
tional comments are in order: (i) Whereas the linear scaling with N of the first moment is 
universal, the scaling of higher moments differ in compact and non compact cases. In partic- 
ular, this scaling implies that while the reduced variance of the FPT is always of order one in 
the non compact case, in the compact case it reads ((T TS ) ~ (Tts) 2 )/ (Tts) 2 oc (R/r) dw ~ d f , 
so very large fluctuations occur for r <C R. (ii) Remarkably, the FPT distribution is entirely 
determined as soon as the mean (Tts) of the FPT is known (as well as its average over the 
starting point (T) T [2"fiH3T)] ) . even if the distribution is not exponential, (iii) For specific 
cases, this mean FPT can be calculated exactly, which provides a fully explicit expression 
of the FPT distribution (as used in fig 2.b, 3. a, 3.b) (iv) In all cases it can be calculated in 
the large volume limit using the recent result [13], leading to the scaling in the geometrical 
parameters given in Eqs(|5])-([6]) . 

We note that our approach covers in particular the important case of subdiffusion [3"T], 
which is characterized by a sublinear dependence of the mean squared displacement with 
time (that is d w > 2). Subdiffusion is widespread in complex crowded environments such as 
biological cells [321 [33] , and might physically originate from a few classes of models based on 
different underlying microscopic mechanisms [34]. Importantly, subdiffusive processes can 
be either compact or non-compact, which will prove below to be the relevant criterion in the 
context of reaction kinetics. The FPT distribution for one class of models for subdiffusion, 
which rely on spatial inhomogeneities [22J as exemplified by diffusion in fractals, is directly 
given by Eqs.([5])-([6]). Another class of models stems from large trapping times, leading to 
the case of infinite (T) T , which we have discarded so far. While the quenched version of this 
type of model becomes quite involved in the case of broad distribution of trapping times 
over the disorder, the FPT distribution in the annealed case - the continuous time random 
walk model (CTRW), which is a standard random walk with random waiting times, drawn 
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from a distribution f(t) [3l | I35H37] - is straightforwardly deduced from Eqs.((5])-(j6]). In this 
case, the Laplace transform of the FPT distribution reads: -Pctrw( s ) — P(f( s ))> where 
P(s) is the generating function (discrete Laplace transform) of the FPT distribution of the 
underlying discrete time random walk, which is determined in Eqs.([5])-([6]). 

These analytical results are validated by Monte Carlo simulations and exact enumeration 
methods applied to various models which illustrate the universality classes defined above. 
These schematic models have been widely used to describe transport in disordered media 
P23E2], for example in the case of exciton trapping on percolation systems [3H] or anomalous 
diffusion in biological cells [39, 40J, as a first step to account for geometrical obstruction and 
binding effects involved in real crowded environments [201 EI]- (i) The non compact and 
marginal cases (see fig [2]) are exemplified by regular diffusion on a 3D and 2D cubic lattices, 
diffusion on a 3D percolation cluster above criticality, and diffusion in disordered systems 
such as the random barrier model (namely a symmetric random walk on a 3D cubic lattice 
with transition rates T drawn from the normalized distribution p(T) oc T~ a ) and the random 
trap model (namely a symmetric random walk on a 3D cubic lattice with frozen waiting times 
T{ at each site drawn from the normalized distribution p(r) oc r~( 1+a >) . (ii) The compact 
case (see fig [3]) is exemplified by diffusion on deterministic fractals such as Sierpinski gasket 
and T-graph and on a critical percolation cluster, as defined in fig 1 of the SI. Figures [2] and 
[3] reveal an excellent quantitative agreement between the asymptotic analytical predictions 
and the numerical simulations, even for systems of small size. We emphasize that despite 
their very different nature, all these models fall into the above defined universality classes. 



Discussion 

We now discuss important implications of these results for reaction kinetics, using the 
ubiquitous example of a target search process involving an immobile target B and searcher 
particles A PQ. When only a small number of A are involved in the reaction, as is the case 

in a microdomain in a biological cell, this reaction has to be described at the single molecule 

t 

level [6] and can be quantified by the survival probability of a particle A, S(t) = 1— P(i'), 

t'=0 

which gives the probability that A has not reacted with B until time t. The quantity S(t) 
depends on the initial position of the molecule A and is explicitly determined using Eqs (§- 
d6]), which indicate that such reactions in confinement obey universal kinetic laws, depending 



7 



on the non compact vs compact nature of the underlying transport process: 

(i) In the non compact case, which corresponds qualitatively to trajectories leaving many 
sites unvisited, as in the case of a 3D medium dilute enough to lead to regular diffusion, 
for any r significantly greater than a typical molecular length, the geometrical factor Uts is 
close to 1, and the dependence on the initial position is lost. We therefore recover a simple 
first order decay of the survival probability which depends on the volume of the confining 
domain only, and not on the initial position of the reactant. In this case of non compact 
exploration, the initial position is not an important parameter of the kinetics (see fig. [4]), 
except in specific cases involving return times, such as recombination reactions . 

(ii) On the contrary, in the compact case, which means that each visited site is on average 
oversampled, geometrical factors dominate. This is typically the case of a crowded medium 
described to a first approximation as a fractal structure where the available space for diffusion 
is restricted. Here, the temporal evolution of S(t) strongly depends on the starting position. 
S(t) drops to small values, indicating that the reaction occurred with high probability, on 
a time scale which depends on the volume, but also critically on the starting position of 
the reactant. This time scale ranges according to EqQ, ^ from r dw (for starting positions 
such that r <C R) to R dw (for starting positions such that r ~ R) , which in practice can 
span several orders of magnitude (see fig. [1]). In this type of "geometry-controlled reactions" 
(not to be confused with "fractal-like reactions" [2]), spatial organization of reactants plays 
a crucial role, which can be quantified by our approach. 

We stress that the decisive criterion leading to geometry controlled kinetics is not the 
subdiffusive vs diffusive nature of the transport process, but its compact vs non compact type. 
We expect this effect to impact a wide class of reactions involving either an inhomogeneous 
initial concentration of reactants, such as a speckeled distribution as experimentally realized 
in |41j . or a small number of particles, such as biochemical reactions in cell subdomains. 
Notably, in the context of gene colocalization, our results give access to the kinetics of 
elementary steps of activation by transcription factors. As an illustrative example, let us 
consider two genes A and B which share a common transcription factor (TF) (for example 
the genes sog and zen of the Drosophila genome which are both targeted by Dorsal [4"2]). 
Experimental results concerning subdiffusive motion of tracer particles in the nucleus [32l 
H31 IB] on the one hand and observations of a fractal organization of the chromatin on the 
other hand (33H3S], provide the following estimates 2 < d w < 3 and df ~ 2.4, and therefore 
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suggest that both compact and non compact exploration cases could occur. Relying on the 
analysis of the survival probability developed previously, we find in the case of compact 
exploration (with for example df = 2.4, and d w = 3 as in |43]) that the typical time needed 
for the TF to reach gene B starting from a gene A colocalized with, i. e. in the vicinity of, 
B (typically tab = r co \ oc < 100 nm, which is the size of a transcription factory [19]) can 
be (rremote/^coioc)^ — 10 3 times faster than for a remote gene A (typically tab = '"remote — 
1 n m, which is the order of magnitude of a nucleus radius). This is in strong contrast with 
the case of non compact exploration (with for example df — 2.4, and d w = 2 [32J) where 
the typical activation time of B starting from A has the same order of magnitude for A 
either colocalized with B or remote. In other words, this suggests that gene colocalization 
is highly favorable for transcription kinetics, but only when it is geometrically-controlled, 
that is in the case of compact exploration, which makes the experimental characterization 
of the nature of transport in the nucleus a major issue. 

To conclude, we calculated analytically the FPT distribution of a diffusing particle to 
a target, and showed that transport processes as various as regular diffusion, anomalous 
diffusion, diffusion in disordered media and in fractals fall into the same universality classes. 
Our results put forward that geometry, and in particular the initial localization of reactants, 
can become a key parameter of reaction kinetics in confinement. In particular, this regime 
of "geometry-controlled kinetics" could be relevant to transcription kinetics and could help 
understand the crucial role of spatial organization of genes. 
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Figure captions : 

Fig 1 : First-passage time distribution (FPT) and geometry controlled kinetics. 

Is (or is not) the initial position of the particle an important parameter of the kinetics ? 
We show quantitatively that in the case of compact exploration (eg for dilute solutions), the 
kinetics turns out to be widely independent of the starting point (Si or S2), whereas in the 
non compact exploration case (eg for crowded environments), the position of the starting 
point strongly influences the search time of the target, leading to "geometry controlled 
kinetics". This result in particular implies that the kinetics of activation of a gene T by a 
transcription factor (TF) can be orders of magnitudes faster if the TF is released from a site 
S = S2 which is colocalized with (i.e. in the vicinity of) T, as compared to the case where 
the TF is released from a remote site S = Si. 

Fig 2 : Universal FPT distribution in the non compact and marginal cases. The 
simulated distribution Gts(9) divided by the weight Uts is plotted against the universal 
theoretical prediction ip(9) Eq(|5]). The collapse of various examples onto a single master 
curve shows the universality of the result, a All non compact and marginal cases (plotted 
independently in b, c, d, e), collapse onto a single exponential master curve, b: Regular 
diffusion on a 3D cubic lattice and 2D square lattice for various rectangular domains (of sizes 
L\ x L 2 x L 3 and Li x L 2 ) and source-target pairs S(x, y, z) and T(x, y, z) whose rectangular 
coordinates are indicated in the legend inset. Here 9 and Uts are calculated using exact 
results for (T TS ) and (T) T given in [21]. c, d, e: examples of disordered systems. Here, 
(Tts) and (T) r are evaluated numerically, c: Diffusion on a 3D percolation cluster above 
criticality embedded in a 30 x 30 x 30 rectangular domain with link probability p = 0.4. d: 
3D random barrier model (see text) embedded in a 30 x 30 x 30 rectangular domain, e: 3D 
random trap model (see text) embedded in various rectangular domains. 

Fig 3 : Universal FPT distributions in the compact case. The simulated distribu- 
tion Gts(8) divided by the weight Uts is plotted against the universal theoretical prediction 
4>(9), Eq([6]). The collapse for different system sizes N shows the universality of the results. 
Examples of deterministic fractals: a Diffusion on a Sierpinski gasket (with a target at the 
apex site) and b on a T graph (with a target at the center). Here, exact expressions are 
used for calculating (Tts) an d (T) T . The insert shows that the scaling function ip weakly 
depends on the dimensions df and d w . c Diffusion on a 3D critical percolation cluster (ran- 
dom fractal) embedded in rectangular domains of sizes (Li x L 2 x L 3 ), as indicated in the 
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inset. Here, (Tts) and (T) T are evaluated numerically, and average over pairs of points is 
performed, to fulfill the scale-invariance hypothesis of the propagator (see text before Eq. 

Fig 4 : Reaction kinetics as quantified by the survival probability S(t), plotted 
for different source-target distances r. The non compact case a is exemplified by a 3D regular 
diffusion and the compact case b by a diffusion on a Sierpinski gasket (with a target at the 
apex). The theoretical prediction for S(t) is obtained from Eqs.([5])-(|6]). Here 9 and Uts 
are calculated using exact results for (Tts) an d (T) T . Quantitatively, the typical reaction 
time t typ can be defined, for example, through the median S(t typ ) = 1/2, indicated by the 
dotted line. In the non compact case, t typ weakly depends on the initial position of the 
reactant, which therefore is not an important parameter of the kinetics. On the contrary, 
in the compact case, t tJP runs over several orders of magnitude depending on the initial 
position which, in turn, controls the kinetics. 
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100x20x20, T(20,5,10), S(70,15,10) 
60x40x20, T(4,4,4), S(56,36,16) 
80x80x10, T(25,25,5), S(55,55,5) 
200x150, T(4,4), S(2,4) 
100x80, T(54,40), S(56,40) 
Theory (Eq. 5) 



random clusters 
and S/T pairs 



Theory (Eq. 5) 



O a=0.8,T(2,2,2), S(28,28,28) 

A a=0.5,T(2,2,2), S(3,3,3) 

□ a=0.2,T(10,10,10), S(20,20,20) 

V a = 0.2,T(2,2,2),S(3,3,3) 

— Theory (Eq. 5) 



O a = 2, 100x20x10, T(2,2,2), S(98,18,8) 
□ a = 1.5, 100x80x10, T(60,60,5), S(20,20,5) 
a= 1.5, 100x80x60, T(3,3,3), S(4,4, 4) 
a = 1 , 1 00x20x20, T(50, 1 0, 1 0), S(40,8,8) 
— Theory (Eq. 5) 
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